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5^ ■ Abstract 

The impact of isolated-photon data from proton-(anti)proton collisions at RHIC, SppS, Tevatron 
and LHC energies, on the parton distribution functions of the proton is studied using a recently 
developed Bayesian reweighting method. The impact on the gluon density of the 35 existing 
isolated-7 measurements is quantified using next-to-leading order (NLO) perturbative QCD cal- 
culations complemented with the NNPDF2.1 parton densities. The NLO predictions are found 
to describe well most of the datasets from 200 GeV up to 7 TeV centre-of-mass energies. The 
isolated-photon spectra recently measured at the LHC are precise enough to constrain the gluon 
distribution and lead to a moderate reduction (up to 20%) of its uncertainties around fractional 
momenta x ~ 0.02. As a particular case, we show that the improved gluon density reduces the 
PDF uncertainty for the Higgs boson production cross section in the gluon-fusion channel by 
more than 20% at the LHC. We conclude that present and future isolated-photon measurements 
constitute an interesting addition to coming global PDF analyses. 
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1 Introduction 

The accurate determination of the parton distribution functions (PDFs) of the proton in a wide 
range of momentum fractions x and energy scales Q [T| is a crucial ingredient for precision phe- 
nomenology at the Large Hadron Collider (LHC) [2l|3]. Among all parton distributions, the gluon 
density g{x,Q'^) is one of the least constrained PDFs since it does not couple directly to the 
photon in deep-inelastic scattering (DIS) measurements of the proton F2 structure function [3]. 
A precise determination of g{x,Q'^) is of paramount importance for the LHC physics program 
both in standard model (SM) as well as in new physics searches. Indeed, gluons drive a sig- 
nificant fraction of the scattering processes at the LHC and gluon-gluon fusion is the dominant 
channel for the production of the SM Higgs boson [5ll6] . top-quark pairs or dijets, to mention a few. 

In global PDF analyses, the gluon density is directly constrained mostly by jet production and 
indirectly constrained by scaling violations of F2{x, Q"^) in DIS and by the momentum sum rule. 
It is thus important to find new collider observables that provide independent information on 
g{x,Q'^) and help improve the accuracy of its determination and reduce its associated uncertain- 
ties. Prompt photon production in hadronic collisions, defined as the production of photons not is- 
suing from the electromagnetic decays of hadrons, appears as an excellent observable to determine 
the gluon PDF, since at leading order (LO) it probes the gluon directly through the quark-gluon 
"Compton" process qg — )• jq [3-l9]. As a matter of fact, up to about 12 years ago, prompt photon 
data were used to constrain g{x,Q'^) in global PDF fits [lOllllj . as well as in independent deter- 
minations of the strong coupling constant |12j . However, a series of measurements carried out at 
centre-of-mass energies y/s w 20-40 GeV by the fixed-target E706 experiment |13H15j showed an 
enhanced 7 cross section compared to next-to-leading-order (NLO) perturbative quantum chro- 
modynamics (pQCD) predictions |16H18| . The theoretical cross section deficit was only partially 
cured by the inclusion of extra soft-gluon resummation contributions |19H23| . Such data-theory 
discrepancies at fixed-target energies, together with the availability of more precise jet measure- 
ments from the Tevatron, lead to abandon the use of inclusive photon data in PDF analyses. 
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The last PDF parametrisation that used the prompt photon data was MRST99 Since 
then, the issue of the compatibihty of NLO pQCD with photon measurements has been discussed 
extensively. For example, the authors of Ref. [21| performed a systematic comparison of NLO 
calculations with most of the available photon results by 2006 (see also Refs. ^5,26j for previous 
studies) and showed that collider isolated photon data was in reasonable agreement with the theo- 
retical predictions. Indeed, as discussed in detail in [2^, by (i) increasing the ^/s from fixed-target 
to collider energies, and by (ii) requiring the photon to be isolated from any hadronic activity 
within a given distance around its direction, one is left with an isolated-photon sample domi- 
nated by energy scales away from non-perturbative effects (such as intrinsic- A;^ broadening |28j). 
and where the uncertain contribution from photons issuing from the collinear fragmentation of 
final-state partons |29] are significantly reduced. Figure [1] shows the relative contribution of the 
three leading diagrams - direct Compton, direct quark-antiquark annihilation, and fragmenta- 
tion - for inclusive (left) and for isolated (right) photon production in proton-antiproton (p-p) at 
Tevatron (top) and in proton-proton {p-p) collisions at the nominal LHC energy (bottom). The 
fragmentation-7 component, including the poorly known gluon-to~photon fragmentation func- 
tion [301 EI] which dominates the low E'^ part of the inclusive photon spectra, is significantly 
reduced after application of isolation cuts and accounts only for less than 10-15% of the cross 
sections at both energies. The right panels of Fig. [1] clearly indicate that the Compton component 
is the dominant contribution to the isolated-7 spectra in wide ranges of the measured transverse 
energy E^, in particular for p-p collisions, thus guaranteeing their increased sensitivity to g{x, Q"^). 
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Figure 1: Relative contributions of the ^(/-Compton, (/^-annihilation and fragmentation subprocesses in 
inclusive (left) and isolated (right) prompt-photon production in p-p (Tevatron at 1.96 TeV, top) and p-p 
(LHC at 14 TeV, bottom) collisions at midrapidity. Results obtained at NLO with jetphox 02] using 
NNPDF2.1 [33], theoretical scales set to I-L = E^, and BFG-II FFs [31] (see Sect.|3|). 
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In view of the various accurate measurements of isolated-photon spectra carried out recently 
at the LHC |34H38j . which have added up about 130 new points to the existing world data, it is a 
timely moment to revisit first if NLO pQCD provides a good description of all available isolated-7 
collider data, and if so, what quantitative constraints these new results impose on the gluon PDF. 
However, the calculation of cross sections at NLO accuracy in hadronic collisions is usually a 
very time-consuming process as it involves the Monte Carlo (MC) generation of a large number 
of event weights for the integration over the hard subprocess phase-space in order to cancel in- 
frared divergences. This fact makes the full-NLO calculation with different PDF choices, as in 
iterative fits of the parton densities, usually prohibitive in terms of the required processing time. 
Computational alternatives exist, though, such as e.g. (i) the ApplGrid framework [39] or the 
FastNLO software [40], both of which make use of lookup tables with cross-section weights which 
can be a posteriori fast-convoluted with any parton density set, and (ii) the NNPDF Bayesian 
reweighting technique |4HI42j which can be applied to any set of MC-based parton densities and 
which relies on the computation of the distribution between the data and each MC replica in 
order to determine the new (reweighted) PDFs. 

In this paper we use the NNPDF reweighting method in order to quantify the impact of col- 
lider isolated-7 data on g{x, Q^). Theoretical NLO pQCD predictions obtained with the jetphox 
program |32j are compared to all available experimental data in a very wide range of centre-of- 
mass (cm.) energies ranging from 200 GeV up to 7 TeV. These comparisons are then used to 
quantitatively determine to which extent isolated-photons can be used to reduce the gluon PDF 
uncertainties. 

The paper is organised as follows. In Sect. [2] the experimental isolated-photon results used 
in the analysis, more than 30 measurements with a total of O (400) data points, are presented. 
Section [3] describes the setup of the calculation based on the the jetphox code to provide NLO 
pQCD predictions and on the NNPDF-reweighting method of the parton densities. The main 
results of this analysis, the data-theory comparisons and the quantitative impact of the photon 
data on the reweighted g{x,Q'^), are discussed in Sect. HI In Sect. [S] the predictions for Higgs 
cross sections in the gluon-fusion production channel using the improved gluon distributions are 
presented. We conclude in Sect. [6] with a summary of the main results and an outline of possible 
future developments regarding the use of isolated-photon data in global PDF fits. 

2 Experimental data 

The world isolated-7 data covers all prompt-photon measurements at collider machines in the last 
30 years. Indeed, for cm. energies above a hundred of GeV, photon isolation is experimentally 
required in order to identify the high-£'^ prompt signal out of the overwhelming background of 
photons from the decays of vr*^ and rj mesons produced in the fragmentation of jets. Table [1] 
summarises the details of the 35 existing experimental measurements. The table lists the char- 
acteristics of the more than 400 isolated-7 data-points measured in p-p and p-p collisions from 
RHIC cm. energies {^/s = 200 GeV) up to the highest energies available so far at the LHC 
{y/s = 7 TeV). The table updates the results collected in [13] by including the latest ATLAS [3S] 
and CMS |38j measurements. For each system we provide the kinematical coverage in E'^ and 
rapidity y-y of the measured photon, the corresponding range of parton fractional momentum x 
covered by the measurement, as well as the applied photon isolation criteria used in each case 
(maximum hadronic energy E/^, or fraction of the photon energy e^, allowed within the isolation 
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cone of radius R = \/ {Uh — 2/7)^ + {4>h — 4>-yY around the photon direction). 



System 


Collab. /Experiment 




\yi\ 


E'^ range 
(GcV) 


X 


Data 


Isolation 




(collider) [Ref.] 


(TcV) 


range 


range 


points 


radius, had. energy 


p-p 


ATLAS (LHC) |34J 


7. 


<0.6 


15-100 


5-10-3-0.05 


8 


R = 0.4, E,, < 5 GeV 


p-p 


ATLAS (LHC) [Ml 


7. 


0.6-1.37 


15-100 


3-10-3-0.1 


8 


R = 0.4, Eft < 5 GeV 


p-p 


ATLAS (LHC) H 


7. 


1.52-1.81 


15-100 


2-10-3-0.1 


8 


R = 0.4, Eft < 5 GeV 


p-p 


ATLAS (LHC) [35] 


7. 


<0.6 


45-400 


5-10-3-0.1 


8 


R = 0.4, Eft < 4 GeV 


p-p 


ATLAS (LHC) ES] 


7. 


0.6-1.37 


45-400 


5-10-3-0.2 


8 


R = 0.4, Eh < 4 GeV 


p-p 


ATLAS (LHC) |35] 


7. 


1.52-1.81 


45-400 


2-10-3-0.3 


8 


R = 0.4, Eft < 4 GeV 


p-p 


ATLAS (LHC) |35] 


7. 


1.81-2.37 


45-400 


2-10-3-0.5 


8 


R = 0.4, Eft < 4 GeV 


p-p 


CMS (LHC) |37] 


7. 


<1.45 


21-300 


5-10-3-0.1 


11 


R = 0.4, Eft < 5 GeV 


p-p 


CMS (LHC) [36] 


7. 


<0.9 


25-400 


5-10-3-0.2 


15 


R = 0.4, Eft < 5 GeV 


p-p 


CMS (LHC) [36] 


7. 


0.9-1.44 


25-400 


2-10-3-0.3 


15 


R = 0.4, Eft < 5 GeV 


p-p 


CMS (LHC) [36] 


7. 


1.57-2.1 


25-400 


10-3-0.4 


15 


R = 0.4, Eft < 5 GeV 


p-p 


CMS (LHC) [361 


7. 


2.1-2.5 


25-400 


10-3-0.5 


15 


R = 0.4, Eft < 5 GeV 


p-p 


CMS (LHC) [38| 


2.76 


<1.45 


20-80 


10-3-0.05 


6 


R = 0.4, Eft < 5 GeV 


p-p 


CDF (Tevatron) [441 


1.96 


<1.0 


30-400 


0.01-0.4 


16 


R = 0.4, Eh < 0.1 


p-p 


DO (Tevatron) [45] 


1.96 


<0.9 


23-300 


0.01-0.3 


17 


R = 0.4, Eft < 2 GeV 


p-p 


CDF (Tevatron) [46] 


1.8 


<0.9 


11-132 


5-10-3-0.2 


17 


R = 0.4, Eft < 4 GeV 


p-p 


CDF (Tevatron) [47] 


1.8 


<0.9 


10-65 


5-10-3-0.1 


17 


R = 0.4, Eft < 1 GeV 


p-p 


CDF (Tevatron) [48] 


1.8 


<0.9 


8-132 


5-10-3-0.2 


16 


R = 0.7, Eft < 2 GeV 


p-p 


DO (Tevatron) [49] 


1.8 


<0.9 


10-140 


5-10-3-0.2 


9 


R = 0.4, Eft < 2 GeV 


P-P 


DO (Tevatron) [49] 


1.8 


1.6-2.5 


10-140 


10-3-0.4 


9 


R = 0.4, Eft < 2 GeV 


P-P 


DO (Tevatron) [50] 


1.8 


<0.9 


9-126 


5-10-3-0.2 


23 


R = 0.4, Eft < 2 GeV 


P-P 


DO (Tevatron) [50] 


1.8 


1.6 - 2.5 


9-126 


10-3-0.4 


23 


= 0.4 Eft < 2 GeV 


P-P 


CDF (Tevatron) [46] 


0.63 


<0.9 


8-38 


0.01-0.2 


7 


R = 0.4, Eft < 4 GeV 


P-P 


DO (Tevatron) [51] 


0.63 


<0.9 


7-50 


0.01-0.3 


7 


R = 0.4, Eft < 2 GeV 


P-P 


DO (Tevatron) [Sl] 


0.63 


1.6-2.5 


7-50 


10-3-0.4 


7 


R = 0.4, Eft < 2 GeV 


P-P 


UAl (SppS) [52] 


0.63 


<0.8 


16-100 


0.03-0.3 


16 


R = 0.7, Eft < 2 GeV 


P-P 


UAl (SppS) [52] 


0.63 


0.8-1.4 


16-70 


0.01-0.4 


10 


R = 0.7, Eft < 2 GeV 


P-P 


UAl (SppS) [52] 


0.63 


1.6-3.0 


16-70 


0.01-0.5 


13 


R = 0.7, Eft < 2 GeV 


P-P 


UA2 (SppS) [53] 


0.63 


<0.76 


14-92 


0.03-0.3 


13 


R = 0.265, £ft < 0.25 


P-P 


UA2 (SppS) [54] 


0.63 


<0.76 


12-83 


0.03-0.3 


14 


R = 0.25, Eft < 0.1 GeV 


P-P 


UA2 (SppS) [51] 


0.63 


1.0-1.8 


12-51 


0.01-0.4 


8 


R = 0.53, Eft < 2 GeV 


P-P 


UAl (SppS) [52] 


0.546 


<0.8 


16-51 


0.03-0.2 


6 


R = 0.7, Eft < 2 GeV 


P-P 


UAl (SppS) [52] 


0.546 


0.8-1.4 


16-46 


0.02-0.4 


5 


R = 0.7, Eft < 2 GeV 


P-P 


UAl (SppS) 'WT 


0.546 


1.6-3.0 


16-38 


0.01-0.5 


5 


R = 0.7, Eft < 2 GeV 


P-P 


PHENIX (RHIC) 55 


0.2 


<0.35 


3-16 


0.03-0.2 


17 


R = 0.5, Eh < 0.1 



Table 1: World systematics of isolated-photon data in p-p and p-p collisions. For each system, we quote 
(i) the experiment, collider and bibliographical reference, (ii) centre-of-mass energy y/s, the measured (iii) 
rapidity \y^\ and (iv) E'^ ranges, (v) the parton fractional momenta x probed, (vi) number of data points, 
and the (v) isolation criteria used. 

Since we have decided not to consider the existing prompt-7 measurements below y/s 65 GeV, 
both at fixed-target energies and at the CERN-ISR collider (see the compilation [25]), as they 
deal with inclusive-^ cross sections with large dependence on poorly known parton-to-photon 
FFs and on non-perturbative corrections due to the smaller scales involved, the lowest-energy 
isolated-7 measurement is that of the PHENIX experiment in p-p at y/s = 200 GeV at RHIC [55]. 
Next in ascending order of collision energies are the oldest measurements by the CERN SppS 
UAl [52] and UA2 |53p54j collaborations at ^/s = 546 and 630 GeV, which amount to a total of 
90 data-points. The Tevatron datasets are available for ^ = 0.63 TeV [MUSI], 1.8 TeV [i6H50] 
and 1.96 TeV [SlIlS], the most precise data being the latter. The total number of data-points is 
21 at 0.63 TeV, 114 at 1.8 TeV, and 33 at 1.96 TeV. Finally, the recent LHC measurements from 
ATLAS and CMS amounting to 133 data points [34ti38] . cover the wider range in rapidity (up to 
|y.y|= 2.5) and have smaller systematical uncertainties (±(7 -20)% depending on and y^) than 
all other previous measurements. For ATLAS we consider both the 880 nb'^ [M] and 36 pb"^ [35] 
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measurements, while for CMS we include the 7 TeV 2.9 pb~^ [37] and 36 pb~^ analyses [36] , as 
well as the latest p-p measurement carried out at 2.76 TeV |38j for reference heavy- ion collisions 
studies. 

The isolated-photon data of Table [Tj are plotted in Fig. [5] as a function of (left panel) and 
as function of = 2E'^l^fs with the cross sections scaled by ^fs"^ (right panel). All spectra in 
Fig. [2] (left) follow clear power-law dependencies from 3 to 400 GeV spanning 9 orders of magni- 
tude in the cross sections. The spectra (Fig. [21 right) coalesce over a single curve when the 
cross sections are normalised by t/s" with exponent n ^ 4.5. Such a behaviour is very close to 
the dependence expected for partonic 2 — )• 2 scattering cross sections in the conformal 

QCD limit, disregarding scaling violations from PDF and running of [56]. A few deviations 
are visible, in particular for measurements that are either away from midrapidity for which the 
assumption = 2E^/y/s does not exactly hold and/or for which the highest photon energies 
are mis-reconstructed (see Sect. 14. ip . The fact that all data satisfy the same scaling with ^/s 
and collapse into a single curve indicates the universality of the underlying partonic production 
mechanism. The observed E^ power-laws and universal -y/s-dependence of the cross sections are 
tell-tale of a perturbative origin for the production of isolated photons and confirm the validity 
of pQCD calculations to study their yields. 
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Figure 2: World systcmatics of isolated-photon spectra measured in p-p and p-p collisions at collider 
energies (Table [Ij as a function of (left) and x.^, (right) where the invariant cross sections have been 
scaled by -^s" with n = 4.5. 

In Fig. [3] we show a scatter plot of the world isolated-photon measurements listed in Table [1] 
together with the DIS, Drell-Yan and jet data (about 3500 points) used in the NNPDF2.1 global 
fits. For each experimental data-point we plot two values of parton fractional momentum x± = 
x^ ■ e^^^/i/i, assuming leading-order partonic kinematics. The plot shows that the photon LHC 
data extend the kinematic coverage in particular over the region x ~ 10~^ — 10^^ at moderately 
large energy scales {E'^'^ ~ 10'^ — 10^ GeV^) not directly covered by the other experimental 
datasets. 
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Figure 3: Kincmatical region in the x ~ plane probed by experimental isolated-7 data at collider 
energies (red circles and triangles) which enter into this analysis (Table [Ij compared to the coverage of 
DIS, Drell-Yan and jet datasets (squares) used in the NNPDF2.1 global fits. 

3 Theoretical setup 

In this section the basic ingredients of the jetphox program used to compute the isolated-photon 
cross sections are discussed, and the NNPDF reweighting technique employed to quantify the 
impact of new data on the proton PDFs is briefly recalled. 

3.1 Isolated-photon cross sections 

Two types of processes contribute at leading order to prompt photon production in p-p and p-p 
collisions: the 'direct' contribution, where the photon is emitted directly from a pointlike coupling 
to the hard parton-parton vertex, and the 'fragmentation' (called also 'anomalous' in the past) 
contribution, in which the photon originates from the collinear fragmentation of a final-state 
parton. Schematically, the differential photon cross section as a function of transverse energy E'^ 
and rapidity can be written as 

da = da^.^ + da^^.^^ = ^ j dxadxf, fa{xa; f4)fb{xb; fJ-l) (1) 

0;b=q,q,9 

/I 
1 Xai X}), Z, fl^ , flp , /J.JJ 
|_ c=q,q,y min 

where fa{xa]fJ^) is the parton distribution function of parton species a inside the incoming pro- 
tons at momentum fraction Xa] d&ab are the parton-parton subprocess differential cross sections; 
and D^/j.{z; fi'^) is the fragmentation function of parton k to a photon carrying a fraction z of 
the parent parton energy, integrated from Zmin = Xrj, coshy^ to 1. The scaled momentum Xrj, is 
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a good representative for the typical x values probed in the PDFs at central rapidities, while at 
forward rapidities one is probing x ~ rc^ exp(— y^). The arbitrary parameters //p., //p and fig are 
respectively the renormalisation, initial-state factorisation, and fragmentation scales which encode 
any residual dependence of the cross sections to higher-order contributions which are missing in 
the calculation. 

The study provided in this article relies on the jetphox calculation of both da,, and da, 

•J ^ dir rrag 

at NLO accuracy [32] in the strong coupling as{Hj^), i.e. all diagrams up to the order 0{aa^) 
are included, in the MS renormalisation scheme. The results of the NLO calculation of da^.^ have 
been known for a long time [57] . The calculation of the NLO corrections to da^^^^ became also 
available later [18 p 32 1 l58]. We note that the distinction between da^^^ and da^^^^ is arbitrary and 
only its sum is physically observable, e.g. bremsstrahlung from a quark leg can be considered as 
"fragmentation" or as "NLO direct" depending on the value of the fragmentation scale considered. 

A few recent works exist that have computed various beyond-NLO corrections to the inclusive 
production of prompt photons: 

• The small-x high-energy corrections j59j were computed in Ref. [60] focusing on the low-E^ 
region of the spectrum where power terms of the type In^(x) are enhanced when the 
scaling variable x becomes small. It was found that such type of next-to-NLO (NNLO) 
corrections are negligible in the kinematical range of the available collider data and, thus, 
ignored in the following. 

• Resummation at next-to-leading-logarithmic (NLL) |19H23| or even next-to-next-to-leading 
(NNLL) [61] accuracy of threshold and recoil contributions due to soft gluon emissions 
which are large close to the phase space boundary when is about half of the cm. energy 
{xrj, —7- 1) have been obtained for da^^^ and da^^.^ . The effect of this resummation is important 
at very large photon corresponding to values Xrj, > 0.2, and provide a much reduced scale 
dependence than the NLO approximation. In the current collider data the typical ranges 
involved {x ~ 0.001-0.2, see Fig. [3]) lead to small threshold corrections which we do not 
consider either. 

• At very large photon transverse energies ~ 1-2 TeV, not yet measured at the LHC, 
one should also consider additional corrections due to electroweak boson exchanges which 
decrease the photon yields by about 10-20% [62] . 

In this paper, we use the jetphox (version 1.3.0) Monte Carlo program |32] to compute the 
NLO pQCD predictions for isolated-photon production. We take into account five active quark 
flavours. The gluon box diagram gg ^ g'j is included in the calculations although its contribu- 
tion to the single inclusive spectrum is found to be just of a few percent. The NNPDF2.1 parton 
densities (with 100 replicas per system) |33j were interfaced to jetphox via the lhapdf (version 
5.8.5) package [63]. The value of as = 0.119 used in the calculations is that provided by the 
NNPDF2.1 parametrisation itself. The renormalisation, factorisation and fragmentation scales 
are all set equal to the photon transverse energy, fi^ = fi^ = = E^. Such fi scales are found 
to result in spectra that agree well with the central values of the experimental data. At LHC 
and Tevatron energies, the sensitivity to changes in the arbitrary theoretical scales was discussed 
in [27] and found to be of about ±(10-15)% in the measured E^ ranges. The parton-to-photon 
fragmentation functions used are the BFG-II ("large gluon") set |31j . The isolated-7 spectrum 
obtained with the alternative BFG-I ("small gluon") set is just a few percent smaller than the 
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one obtained with BFG-II in the lowest range [27] since a significant fraction (up to 80%, see 
Fig. [1]) of the fragmentation photons are removed by the application of isolation cuts. 

We run JETPHOX for all the systems in Table [1] within the {E^,y-y) kinematics ranges of every 
measurement. The MC photon isolation criteria are also matched as closely as possible to each 
one of the experimental cuts. We histogram the partonic configurations generated in S^-bins of 
the same size as for the experimental data in order to be able to compute the bin-by-bin and 
avoid potential problems with steeply- falling spectra such as those discussed in Ref. [6l]. Running 
for each one of the 100 NNPDF2.1 replicas with acceptable MC statistics (statistical uncertainty 
below 1%) for all bins is a very CPU-demanding task. The jetphox MC production for 100 
PDF replicas takes about 7 (resp. 10) hours per each million direct (resp. fragmentation) events 
in a O (2 GHz) CPU core, which means that about one week of computer-time is needed to obtain 
the 100 spectra for each one of the 35 systems considered. 



3.2 PDF reweighting 

In order to quantify the impact of the collider isolated-photon data on the proton PDFs we use 
the Bayesian reweighting method described in Refs. |4HI42j. This technique can be applied to 
any PDF-set that estimates PDF uncertainties based on a MC method, such as the NNPDF 
family [MlllllESHZQ] (but also MSTW [7l] and HERAPDF [72] have produced MC sets). The 
method allows one to straightforwardly determine the impact of a new dataset on PDFs by means 
of Bayesian inference. The only ingredient needed are the data-theory "goodness-of-fit" x| for 
each k-th MC replica (where k = l,...,iVrep runs over the full set of replicas: A'rcp = 100 or 1000 in 
the NNPDF2.1 set). This technique reduces the problem of the slowness of jetphox (or in general 
any NLO hadronic computation) which forbids its direct use within a PDF-fitting program, since 
the theoretical predictions need to be computed only once for each PDF replica. 

The Xk each replica is defined as follows: 



xl 



1=1 



where af"^^'^^^ is the NLO theoretical prediction for the isolated-photon cross section, Eq. ([T|), 
obtained with the fk PDF replica, a^^^^^ is the corresponding experimental measurement and 
Atot accounts for the experimental uncertainties. A'^jjat is the number of data points of each par- 
ticular measurement. We note that the experimental covariance matrix is not available for any 
of the isolated-photon datasets and thus we are forced to add in quadrature the statistical and 
systematic uncertainties into Atot- Likewise, absolute normalisation uncertainties (e.g. from the 
integrated luminosity) should be in principle included in the experimental covariance matrix us- 
ing the to method (as consistently done within the NNPDF analysis) [73], but they are typically 
smaller than other experimental uncertainties and also added in quadrature here. 

Once the xt ^or each replica has been computed, the new weight of the replica is given by 

(^2)i(n-l)g-ixl 
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The weights Wk, when divided by the number of MC rephcas of the prior PDF set (A'^rcp)i are then 
simply the probabihties of the rephcas fk, given the xt added experimental results. If 

the new data constrains the PDFs, reweighting will be less efficient that refitting because of the 
discarded replicas with low weight. One can quantify this efficiency loss by using the Shannon 
entropy to compute the effective number of replicas left after reweighting: 

{JVrcp "\ 
jt^Y^Wk ln(iVrcp/wfc) > . (4) 

Clearly < N^s < N^cp and the reweighted fit has the same accuracy as a refit with N^s replicas. 
This effective number of replicas A'gff is a useful estimator of the impact of each individual dataset. 
The smaller A'^eff is, the more the new dataset constrains the PDFs. 

It is often the situation that experimental uncertainties are under or over estimated. It is then 
possible to rescale the uncertainties of the data by a factor a, and then use inverse probability to 
calculate the probability density for the rescaling parameter a 

Nrep 

k=l 

Here Wfc(a) are the weights Eq. ^ evaluated by replacing xl with xl/a^, and are thus propor- 
tional to the probability of fk given the new data with rescaled errors. Averaging Wk{a) in the 
reweighted fit thus gives the probability density for a. If this probability density peaks close 
to unity, the new data are consistent with the pre-reweighting data, while if it peaks far above 
(below) one then it is likely that the errors in the data have been under (over) estimated. The 
distribution of the rescaling variable a, normalised to unity, is used in the next section to inves- 
tigate if the bad agreement between NLO pQCD and a few datasets could be due to possibly 
underestimated experimental uncertainties. 



4 Results 

The comparison between NLO pQCD and the experimental data for all systems listed in Table [T] 
is presented in the next subsection. With this information we compute for each case the data- 
theory Xk distributions Eq. ([2]), the weights Eq. ([3]), the effective number of replicas Eq. (j4]), 
the sensitivity to individual PDF flavours, and finally the associated constraints on the gluon 
distribution. 



4.1 Comparison between data and NLO pQCD 

Using the jetphox program with the setup discussed in Sect. 13.11 we have computed the E^- 
differential cross sections using the A'rcp = 100 replicas of the NNPDF2.1 NLO parton distri- 
butions for all the systems listed in Table [U The ratios between data and theory are shown in 
Figs. IIHHI where the (yellow) band corresponds to the distribution of the predictions obtained 
which each one of the 100 replicas, while the outer error-bars cover the sum in quadrature of the 
statistical and systematic uncertainties of the measurement. 

Figure [H shows the data/NLO ratios for the isolated-photons measured at the lowest cm. 
energies: ^/s = 200 GeV (RHIC) and 546 GeV (SppS). The PHENIX results agree well with NLO 
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pQCD within the quite large experimental uncertainties. For UAl and UA2, the data/theory 
ratio is around unity below w 25 GeV, but the highest transverse energy points are clearly 
overestimated by NLO in two cases, likely due to an experimental mis-reconstruction of the pho- 
ton energy at the end of the spectra measured at the time (see below). 
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Figure 4: Ratio of isolated-photon data and NLO pQCD predictions for p-p collisions at y/s = 200 GcV 
(PHENIX) and p-p collisions at ^/s = 546 GeV (UAl) and various photon rapidities. The (yellow) band 
indicates the range of predictions for each one of the 100 NNPDF2.1 replicas, and the bars show the total 
experimental uncertainty. 



The comparison of the theoretical predictions to the UAl, UA2, CDF and DO data at 630 GeV 
cm. energies and various photon rapidities are shown in Fig. [5j With a few exceptions, a general 
trend appears indicating that the measured cross-sections are higher (respectively lower) than the 
NLO calculations at the low (resp. high) end of the photon energy spectrum. The disagreement 
is larger for increasingly forward rapidities. Yet, given the relatively large experimental uncer- 
tainties, theory and data are in general consistent and show not very far from one in most of 
the cases. 



The Tevatron Run-I measurements at ^/s = 1.8 TeV are compared to the NLO predictions 
in Fig. [6l The two oldest measurements |481l50j (not plotted here) are systematically below not 
only the pQCD predictions but also the spectra obtained by the same experiments in a basically 
identical kinematic range a few years later |461I47^ H9] . By discarding those superseded datasets 
from further analysis, we find a generally good consistency of the 1.8- TeV data with NLO pQCD. 

The comparison of NLO to the Tevatron CDF and DO Run-II results at ^/s = 1.96 TeV is 
shown in Fig. [71 While at the level of the data-theory agreement is good, in both cases 
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Figure 5: Ratio of isolated-photon data (UAl and UA2 in the 3 top and middle plots, and CDF and DO 
in the 3 bottom panels) and NLO pQCD predictions for p-p eollisions at ^/s = 630 GeV at various photon 
rapidities. The (yellow) band indicates the range of predictions for each one of the 100 NNPDF2.1 replicas, 
and the bars show the total experimental uncertainty. 
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Figure 6: Ratio of isolated-photon data (CDF and DO) and NLO pQCD predictions for p-p collisions at 
-y/s = 1.8 TeV at various photon rapidities. The (yellow) band indicates the range of predictions for each 
one of the 100 NNPDF2.1 replicas, and the bars show the total experimental uncertainty. 



the spectral shape seems somewhat different at small photon E^j, where the data rise steeper as 
compared to the theory. The origin of this different shape is still not well understood, although 
Ref. [23] argues that the discrepancy decreases with a suitable scale choice. 
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Figure 7: Ratio of isolated-photon data (CDF and DO) and NLO pQCD predictions for p-p collisions at 
i/s = 1.96 TeV at central rapidities. The (yellow) band indicates the range of predictions for each one of 
the 100 NNPDF2.1 replicas, and the bars show the total experimental uncertainty. 

Data-theory comparisons for ATLAS and CMS measurements are shown in Fig. [8] and [9] re- 
spectively. We observe a very good agreement for all rapidity ranges, except maybe in a few of the 
lowest-£'^ bins where the central value of the data points tends to undershoot a bit the theoretical 
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predictions. 
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Figure 8: Ratio of ATLAS isolated-photon datasets (880 nb~\ top, and 36 pb~\ bottom) and NLO pQCD 
predictions in p-p collisions at ^/s — 7 TeV at various photon rapidities. The (yellow) band indicates the 
range of predictions for each one of the 100 NNPDF2.1 replicas, and the bars show the total experimental 
uncertainty. 
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Figure 9: Ratio of CMS isolated-pliotoii data and NLO pQCD predictions in p-p collisions at 
Vs = 2.76 TeV (top left) and = 7 TeV (2.9 nb~\ top-middle, and 36 pb-\ top-right and bottom, 
datasets) at various photon rapidities. The (yellow) band indicates the range of predictions for each one 
of the 100 NNPDF2.1 replicas, and the bars show the total experimental uncertainty. 

Figure [TOl shows the overall comparison of the full collider photon dataset considered in this 
analysis with the NLO predictions obtained using the central NNPDF2.1 replicas, as a function 
of Xrj,. Most of the data/theory ratios are around unity indicating an overall good agreement 
between NLO pQCD and the experimental measurements. No systematic deviation is observed 
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in a wide range of x^j, corresponding to a wide kinematical coverage of photon , rapidity and 
collision energy. This is to be contrasted with similar systematics studies |24j that indicate clear 
data-NLO deviations for the inclusive-7 production in the E706 results at fixed-target energies. 




Figure 10: Summary plot for the data/theory ratios for all collider isolated-photon data considered in this 
analysis, as a function of = 2E^/^/s. For each system, the NLO prediction used is the one obtained with 
the central PDF replica of the NNPDF2.1 set. The error bars indicate the total experimental uncertainty. 

In Table [1] we quote the average over all replicas between each one of the datasets and the 
NLO calculations. As one can see, for a large majority of cases the agreement is quite good, while 
in a few cases the obtained is rather poor (x^ ^ !)■ The total initial of all the systems 
considered is x'^ = 1-3, while after reweighting it decreases to Xrw = 1-1 for the whole dataset. 
This global result confirms at the quantitative level that there is a good agreement between NLO 
pQCD and the experimental results measured at all collider energies. 

Looking in more detail, it is worth noticing that for those few systems which are not well 
reproduced by NLO pQCD, there are always other measurement^ covering similar {^,E^,y^) 
domains with a good x^- This fact indicates that in such cases the problem is not likely re- 
lated to the theoretical prediction but of experimental origin. In other words, the issue is not 
data-theory compatibility but rather and inconsistency problem between measurements covering 
the same kinematics. Note that in the cases where the x"^ is poor, the reweighted Xrw only 
slightly better, which confirms that the data are not consistent with the theory even after refitting. 
To quantify better this effect, we can make use of the probability distribution for the rescaling 
variable a discussed in Sect. 13.21 The last column of Table [2] lists the mean value of the V{a) 
distributiorll, Eq. ([5]), for all systems. In Fig. [TT]we show the V{a) distribution for those systems 

^The only exception to this are two UAl measurements at 546 GeV (Fig. 3)) for which the large is just driven 
by a single outlier data-point at the highest measured. 

^It is easy to check that if the underlying distribution is a distribution, the mean value of 'P(a) is given by 
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Table 2: Summary of the x^-analysis between NLO pQCD and the world isolated-7 data. For each 
system we list the initial data-theory (6th column), the x?w (7th column) obtained after including 
each corresponding dataset via PDF reweighting, and (a) (last column), the mean of the associated 7^(0!) 
distribution. 
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in Table [2] with initial ^ 3. In all cases V{a) peaks at a value of a above ~2. This fact points 
to an underestimation of the experimental uncertainties in the corresponding measurements. 




Figure 11: Distribution of the a rescaling variable, Eq. (O, for those datasets that show poor agreement 
with NLO pQCD. From top-to-bottom and left-to-right, we show two UAl-546 GeV datasets (x^ = 6.1 and 
7.4 respectively), one UAl (x^ = 5.9) and UA2 (x^ = 3.0) datasets at 630 GeV, and finally two Tevatron 
Run-I CDF and DO sets (x^ = 3.6 and 3.1 respectively). 

As a comparison, in Fig. [T2]we show the 'P(a) distributions for the rapidity bins of the 36 pb~^ 
results from ATLAS and CMS. In all cases V{a) peaks close to one, confirming the consistency 
of these datasets with NLO pQCD and the proper estimation of their associated experimental 
errors. 




Figure 12: Distribution of the a rescaling variable, Eq. ([5]), for each one of the rapidity bins of the 36 pb 
datasets from ATLAS (top plots) and CMS (bottom plots). 



(a) = 1-1- l/2A'^dat with A'dat the number of data points of that particular system. 
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4.2 Sensitivity of isolated-7 to the proton PDFs 



The dependence of isolated-7 production at different kinematics ranges on the individual flavour of 
the underlying parton densities can be quantified by computing the correlation coefficient between 
each one of the light-quark and gluon distributions and the NLO cross sections as discussed in [33j . 
Correlation profiles are shown in Fig. [13] for four representative values of the photon energy and 
rapidity typical of the LHC measurements. Isolated photons at central rapidities (top panels) 
have a dominant sensitivity to g{x,Q'^) in a range between x = 0.01 and x = 0.2 for increasing 
photon transverse energy. At forward rapidities (bottom panels) isolated-7 are sensitive to the 
gluon densities at x values around 10^"^ for low E^, but the PDF sensitivity shifts to the valence 
it-quarks as the photon energy increases. Being able to carry out measurements at even more 
forward rapidities, e.g. down to ~ 5 accessible in LHCb [71] and/or at low-E^ accessible via 
photon-conversions in ALICE [75], would be very interesting in order to probe x values below 
10~^ where g{x,Q^) is only weakly constrained [76] by HERA F2, Fl, F2.C structure functions 
and by the momentum sum-rule. 
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Figure 13: Correlation between the isolated-7 cross section and various flavours of the NNPDF2.1 parton 
densities for different {E^,y^) kinematical ranges at the LHC. Top (bottom) plots show the correlation at 
central (forward) rapidities. Left (right) plots show the correlation for typical low (high) values. 

In Fig. [13] we plot the sensitivity of the isolated-7 data to the proton gluon at a fixed value 
of the photon transverse energy (E^ = 27 GeV) for various collider energies at central (left) and 
forward (right) rapidities. At midrapidity, the peak of the correlations shifts towards larger x 
values for decreasing ^/s, up to about x ~ 0.1 at O (500 GeV). At forward rapidities the LHC 
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data are more sensitive to the gluon density than measurements at SppS and Tevatron because 
the former involve p-p colHsions (dominated by g^f-Compton processes) rather than p-p colhsions 
(with gg-annihilation playing a larger role, Fig. [T]), and because the {E'^,y^) phase-space covered 
by ATLAS and CMS is much larger than that at lower cm. energies. 
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Figure 14: Correlations between the gluon PDF and the isolated-7 cross section at a fixed = 27 GeV 
in p-p and p-p collisions at central (left) and forward (right) rapidities for different collider energies. 



4.3 Impact of isolated-photon data on the gluon 

In order to quantify the constraints that the various photon datasets impose on the gluon PDF, 
the results of the x^-analyses listed in Table [2] are combined into seven groups according to each 
collider energy : 200, 546, 630 GeV and 1.8, 1.96, 2.76 and 7 TeV. In the grouping, we discard the 
six datasets with ^ 3 which have underestimated experimental uncertainties (Fig. [TT]) and only 
result in a loss of accuracy of the reweighting method without any impact on the PDFs whatsoever. 

The effective number of replicas A'eff after reweighting, Eq. Q, for each collider energy is 
summarised in Table [3l As expected, the more constraining datasets, i.e. those with smallest 
N^Q, are those from the LHC at 7 TeV. The measurement at 2.76 TeV has only 5 data-points 
fully consistent with NLO with basically no impact on the PDFs. Of the rest, only the Tevatron 
Run-II and the 630 GeV datasets seem to have also some, albeit small, reduction of A'^cff- 

y/s (TeV) (12 0.546 0.630 L8 L96 2?76 ~ 
A^eff 99.6 99 95 99.8 96 96 87 

Table 3: Effective number of replicas Ngg, Eq. (jlj, left for each one of the collider energies considered, after 
inclusion of the isolated-7 data. The starting number of NNPDF replicas for this analysis is A^iop = 100. 

In Fig. [15] we show the 100 replica weights, Eq. ([3]), for three representative collider ener- 
gies: PHENIX at y/s = 200 GeV, UAl/UA2/Tevatron at ^ = 630 GeV, and the LHC data at 
^/s = 7 TeV. For the least constraining dataset (the PHENIX data with N^g closer to A^ep) all 
replicas have essentially the same weight. For the other two datasets the distributions become 
broader, specially for LHC at 7 TeV, showing that the impact of these measurements on the 
PDFs is larger - some replicas are preferred (larger weight) than others (smaller weight) when 
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confronted to the photon results. 
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Figure 15: Distribution of NNPDF replicas weights, Eq. ([3|), for three representative collider energies: 
PHENIX data at = 200 GcV (left), = 630 GeV data (center), and LHC results at = 7 TcV 
(right). 



The direct quantification of the impact on the g{x,Q'^) distribution is shown in Figs. [TBHT51 
where the NNPDF2.1 NLO gluon is shown before and after including the isolated-photon data 
with reweighting. In Figs. [TBHTTl the PDFs are evaluated at a typical LHC scale oi Q = 100 GeV. 
We show both the relative improvement (left panels, ratio over the original gluon PDF) as well 
as the reduction in the absolute PDF errors (right panels). Only the LHC data lead to a sig- 
nificant PDF uncertainty reduction, of up to 20 percent, localised at medium x ~ 0.002 to 0.05 
(Fig. [T8]l . The Tevatron Run-H and the 630-GeV measurements bring rather small improvements 
around x ~ 0.01-0.02, while the other datasets have negligible impact on the gluon central value 
and associated uncertainties. This is consistent with the expectation that the effects should be 
maximal in phase-space regions where the photon production cross section depends more strongly 
on the gluon distribution (Fig. [T4|) . In the cases where the PDF uncertainties are reduced, the 
central value of g{x,Q'^) is essentially unaffected (Figs. [T6HT71 left) indicating that the large-x 
gluon determined from the Tevatron jet data is consistent with the large-x gluon constrained by 
the LHC photon results. This is an important cross-check of the validity of pQCD factorisation 
and of the PDF universality using cross sections measured at hadronic colliders. The effect of the 
LHC data on the gluon is similar at lower scales, where the impact of photon data is shifted to 
somewhat larger values of x as dictated by DGLAP evolution, as exemplified in Fig. [181 

The relative reduction of the uncertainties in the gluon distribution is better illustrated in 
Fig. [El For three scales Q = 3.16, 10, and 100 GeV, the maximum g[x,Q'^) uncertainty reduc- 
tion is 10-20% in the region x w 0.01-0.05. For the lower scales, in some small regions of x, the 
PDF errors can be slightly increased when the photon data are included. This is just a fluctuation 
effect, due to the limited statistics of the reweighting procedure, that is quickly washed out by 
DGLAP evolution at higher scales. 

As a last check, we have confirmed that the constraints on the quark PDFs from the current 
isolated-7 data are essentially negligible. This is not unexpected since light-quarks distributions 
are known more accurately than the gluon PDF in the kinematical region relevant for photon 
production (Fig. [T3l) . However, the impact of the isolated-photon measurements would certainly 
be more important in the so-called "collider only" fits j70j . where the use of data from colliders 
alone leads to larger PDF uncertainties on the quark sector (constrained mostly by fixed-target 
data in pre-LHC PDF sets), since in this case the q, q densities could be also directly constrained. 
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Figure 16: Comparison between the NNPDF2.1 NLO gluon before (green solid band) and after (dashed 
blue area) inclusion of the isolatcd-7 data from (top to bottom): LHC-7 TeV, Tevatron Run-II at 1.96 TeV, 
and Run-II at 1.8 TeV. The left plots show the ratio between the original and the new Q^) while the 
right panels indicates the reduction of absolute g{x^ Q^) uncertainties thanks to the photon data. PDFs 
are valued at Q = 100 GeV, a typical LHC scale. 
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Figure 17: Same as Fig. [16] for the remaining collider data (top to bottom); UAl/UA2/Tevatron at 
630 GeV, UAl-546 GeV, and PHENlX-200 GeV. 
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Figure 18: Comparison between the NNPDF2.1 NLO gluon before (green solid band) and after (dashed 
blue area) inclusion of the LHC-7 TeV isolated-7 data, with PDFs evaluated (left to right) at Q = 3.16, 
10, and 100 GeV. 



LHC 7 TeV isolated-y data 

^ 1-3r 

cvi - NNPDF2.1 NLO + lsoPhotons(Q^=10GeV^) 

u. 

g 1.2- NNPDF2.1 NLO + lsoPhotons(Q^=100GeV^) 

^ : NNPDF2.1 NLO + lsoPhotons(Q^= 10000 GeV^) 

o 1.1- 




'■^ 0.0002 0.001 0.010.02 0.1 0.2 



X 

Figure 19: Relative reduction of the NNPDF2.1 NLO gluon distribution uncertainty at scales Q = 3.16, 
10, and 100 GeV, after inclusion of the LHC isolated-photon data via reweighting. 
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5 Predictions of the gluon-fusion Higgs boson production cross 
section 

The dominant production channel for the SM Higgs boson at the LHC is gluon fusion OE] and 
thus theoretical predictions depend strongly on the gluon PDF choice as well as on the associ- 
ated value of as- Recently there has been an intense discussion in the literature |77fl79] about 
the theoretical uncertainties to assign to the cr(g g — )■ H) cross section at hadron colliders. This 
has been motivated by the fact that predictions from non-global PDF sets lead to very different 
results than global PDF fits which are in reasonable agreement among each other [2] . This has of 
course implications for the current Tevatron and LHC Higgs exclusion limits |8UH82j . Dedicated 
studies |83H85j have shown that, within a global fit, the gluon PDF and the strong coupling are 
stabilised by the inclusive jet data. It is thus of outmost importance to find additional observables 
sensitive to g{x,Q'^) to further improve the predictions for gluon-gluon Higgs cross sections, and 
in this context isolated-photons appear as a promising candidate. 

As shown in the previous section, the LHC isolated-photon data leads to a reduction of the 
g{x,Q'^) uncertainties at around x ~ 0.02. Since this kinematical region is relevant for many 
important gluon-driven processes at the LHC, we study how the PDF uncertainties involved in 
the determination of the cross sections for various processes, from Higgs boson to top-pair produc- 
tion [86], can be potentially reduced thanks to the inclusion of isolated-7 data. We have decided 
to consider only 7-TeV photon results since the constraints from measurements at the other cm. 
energies are much milder and, in addition, LHC is the only collider for which no dataset needs 
to be discarded. We show also that the central cross section predictions for Higgs production, 
mostly driven by the inclusive jet data in global PDF fits, remain stable when isolated-photon 
data are included. 



Process / Cross section gg g(120) gg g(160) gg H(200) gg g(500) 

NNPDF2.1 11640 ± 181 fb 6052 ± 103 fb 3494 ± 66 fb 219.3 ± 8.3 fb 

NNPDF2.1 + LHC IsoPhotons 11701 ± 140 fb 6073 ± 86 fb 3504 ± 56 fb 218.4 ± 7.6 fb 



Process / Cross section ti_ ttH{120) W H{12Q) Z H(12Q) 

NNPDF2.1 162 ± 51 pb 114 ± 5 fb 447 ± 9 fb 364 ± 6 fb 

NNPDF2.1 + LHC IsoPhotons 162 ± 47 pb 113 ± 4 fb 448 ± 9 fb 365 ± 6 fb 

Table 4: NLO cross sections in p-p collisions at 7 TeV obtained with mcfm using the NNPDF2.1 set, 
before and after including the LHC isolated-photon data. Top: Higgs production in gluon-gluon fusion 
for different values of A/h- Bottom: QCD top-pair production, and Higgs production (Mh = 120 GeV) 
associated with top-pairs and electroweak bosons. 

The production cross sections for the following processes have been computed at NLCH with 
the MCFM code [87ti89] for p-p collisions at \/s= =7 TeV: (i) Higgs boson in gluon fusion {gg — )• H) 
for masses Mh = 120, 160, 200, 300 and 500 GeV, Higgs boson in association with (ii) top quark 
pairs (ttH), and (iii) W and Z {WH, ZH); and (iv) top quark pair production (tt). All these 
processes are, either as signal or as background, relevant for Higgs searches at the LHC. The 
results for the cross sections obtained with the NNPDF2.1 NLO set, before and after including 
the LHC isolated photon data, are summarised in Tabled! The results for Higgs boson production 

^The NNLO/NLO K-factor, when available, is know to depend mildly on the choice of PDF. Therefore, for the 
study of interest here (impact of the PDF error reduction on LHC cross sections) NLO calculations are sufficient. 
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in gluon-gluon fusion are also shown in Fig. [20l While the most updated searches exclude the SM 
Higgs boson with masses outside the 115-127 GeV range |8HI82]. we nevertheless show a wider 
mass range to illustrate the impact of the photon data. 
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Figure 20: Ratio of Higgs production gluon-fusion cross sections with NNPDF2.1 NLO PDFs before and 
after including the LHC isolated-photon data. 



As one can see from Table H] and Fig. [20l the inclusion of isolated-photon results in the NNPDF 
analysis leads to reduced PDF uncertainties in the gluon-gluon Higgs production cross section, 
which is maximal for low Higgs masses (in the region not yet excluded by the ATLAS and CMS 
limits) and can be as large as 20%. For other processes the improvements are much more mod- 
est. Indeed, processes with a pair of top-quarks are produced with much larger virtualities than 
a (low-mass) Higgs and have reduced gluon uncertainties. The same holds for processes with 
associated electroweak boson production which are dominated by quark PDFs. In all cases, the 
central prediction is in good agreement with the reference NNPDF2.1 results, and thus consistent 
with the information obtained from the inclusive jet data in a global fit. All in all, these results 
demonstrate that isolated-photon data can be useful to reduce the theoretical (PDF) uncertainties 
in the cross section for the many gluon-induced processes at the LHC. 
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6 Summary and outlook 



Using up-to-date NLO pQCD theoretical calculations for isolated-photon production, as imple- 
mented in the jetphox program combined with the NNPDF2.1 parton densities and its associated 
PDF reweighting technique, we have quantified the impact on the gluon density of all the existing 
isolated-7 data measured in p-p and p-p collisions at collider energies. The main results of this 
work can be summarised as follows: 

• NLO pQCD provides generally a good description of the isolated-photon measurements from 
200 GeV up to 7 TeV in a wide kinematic range of photon transverse energies and rapidities. 
A few outliers exist which can be identified as arising from problems in experimental datasets 
with underestimated systematic uncertainties. 

• The ATLAS and CMS measurements at the LHC, amounting to more than one hundred 
data-points, are consistent with each other and with pQCD and precise enough to mod- 
erately constrain the gluon PDF in the region x ~ 0.002 — 0.05 in a large range of scales 

= 10-10^ GeV^. The central value of the gluon distribution is unmodified but its un- 
certainty is reduced by up to 20% around x « 0.02. 

• New PDFs including the LHC photon data lead to improved predictions for low mass Higgs 
production in the gluon-fusion channel, with central values of the theoretical cross section 
unmodified but with associated PDF uncertainties decreased by as much as 20%. 

These results confirm that there is no reason not to include isolated-photon data measured at 
LHC energies into global-fit PDF analyses, on a similar footing as the current inclusive jet pro- 
duction data. The main difficulty to systematically carry out such a program is the slowness of 
the NLO pQCD codes to compute the corresponding theoretical predictions, and thus an eventual 
fast interface d la ApplGrid or FastNLO would be very advantageous. 

Various possible future developments of this study are outlined next. First, from the exper- 
imental side, measurements of the inclusive single isolated-photon spectra at high-i?^ exploiting 
the much larger luminosity collected in 2011 by the ATLAS and CMS collaborations at 7 TeV, as 
well as data at 8 TeV in 2012 and eventually at 14 TeV would be very beneficial. The addition 
of a few hundred new data-points over a large kinematic range, possibly including the enhanced 
E^-y^ range covered by ALICE and LHCb, will certainly lead to even stronger constraints on 
the gluon density. The measurement of ratios of photon cross sections at 8 (or 14) to 7 TeV, 
where part of the experimental uncertainties cancel out, and/or of the double-differential 7-jet 
spectra (with constrained kinematics by the concurrent measurement of the photon and jet) can 
also provide new important constraints on g{x, Q^). Second, at the lowest end of collider energies, 
the possibility of measuring isolated-7 production in the full range of the RHIC cm. energies 
(y^ ?a 20-500 GeV) would be also very useful to clarify for once the long-standing disagreement 
between fixed-target data (as well as the oldest ISR and a few of the SppS results) and NLO calcu- 
lations. In all cases, the use of a smooth-cone prescription for the photon isolation [90] which can 
effectively remove any remaining fragmentation-photon component would be helpful to further 
eliminate theoretical uncertainties. Third, it would be also important that future photon results 
provide the full covariance matrix, just as recent LHC measurements of electroweak boson [91] 
and inclusive jet [92] production have done, as this would improve the quantitative treatment of 
these datasets in PDF analyses. 
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On the theoretical side, it would be useful to carry out a more quantitative study of the 
possible impact of scale variations on the photon cross sections. It is conceivable that part of 
the impact of the photon data is washed out once higher-order uncertainties are accounted for. 
This situation is common to all hadronic production data from colliders (notably, inclusive jets) 
included in global fits. One possibility to reduce this problem would be to include threshold re- 
summation corrections. In any case, such a scale variation analysis should be consistently carried 
out together with all the other datasets. Once photon data are included in PDF analyses, it could 
be also important to accurately and consistently determine the strong coupling in the framework 
of a PDF fit, just as inclusive jet production is now instrumental in that regard |83 P85l[93] as 
compared to DIS-only fits. This is so because the gluon is loosely constrained in DIS-only fits, 
where it has a runaway direction. 

As a conclusion, we have shown that the available isolated-photon data provides constraints 
on the gluon PDFs and thus on many relevant LHC processes, most importantly Higgs production 
in gluon-gluon fusion. Given that even more precise data as well as theoretical improvements will 
be available in the next future, we see no objection why isolated-photon data should not become 
integral part of future global QCD analyses. 
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